Breakdown of Diffusion in Dynamics of Extended Waves in Mesoscopic Media 
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We report the observation of nonexponential decay of pulsed microwave transmission through 
quasi-one-dimensional random dielectric media that signals the breakdown of the diffusion model 
of transport for temporally coherent extended waves. The decay rate of transmission falls nearly 
linearly in time due to a nearly gaussian distribution of the coupling strengths of quasi-normal 
electromagnetic modes to free space at the sample surfaces. The peak and width of this distribution 
scale as L -2 05 and L -1 ' 81 , respectively. 
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The diffusion model is widely applied to electronic, 
neutron, and thermal conduction as well as to acoustic 
and electromagnetic propagation in multiply scattering 
media. The model is used not only when the phase of the 
wave is scrambled by inelastic scattering, but also when 
the wave is temporally coherent in mesoscopic samples 
P, 0] • Although wave interference leads to large inten- 
sity fluctuations within a particular sample, the ensemble 
average of the flux reaching a point is generally assumed 
to be the incoherent sum of contributions of randomly 
phased, sinuating Fcynman paths reaching that point. In 
this model, the average intensity varies smoothly in space 
and time and is governed by the diffusion equation. 

Diffusion has been taken as the counterpoint to wave 
localization [1- [J 0] On one side are sharply defined 
localized modes with the average level spacing Av ex- 
ceeding the typical level width bv and on the other are 
diffusing waves, for which bv > Av with bv ~ D/L 2 , 
where D is the diffusion coefficient and L is the sample 
thickness [jj. However, the mode picture is inescapably 
a wave picture and at variance with particle diffusion in 
a number of respects. First, the diffusion equation is of 
first order, whereas the wave equation is of second order 
in time. The time evolution of the wave at any instant 
should therefore depend not only upon the spatial distri- 
bution of the particle density or intensity at that instant, 
as it does in the particle diffusion picture, but also on the 
previous history of the wave. Second, the particle picture 
represents the intensity as a discrete sum over diffusion 
modes, whereas the wave picture describes the field as a 
superposition of quasi-normal modes with a continuum of 
decay rates in a random ensemble. The decay rates of the 
diffusion modes are given by l/r„ = n 2 n 2 D/(L + 2z ) 2 , 
where n is a positive integer and zq is the boundary ex- 
trapolation length. After a time ti, the intensity distri- 
bution settles into the lowest diffusion mode and decays 
at a constant rate, 1/t±. In contrast, the decay rates of 
quasi-normal modes in a random ensemble should be a 
continuum. As time progresses, long-lived quasimodcs 
would contribute more substantially and the rate of flow 
out of the sample would slow down continuously. Nonex- 



ponential decay has been observed in acoustic scattering 
in reverberant rooms 0] and solid blocks Q as well as 
in microwave scattering in cavities whose underlying ray 
dynamics is chaotic . Similarly, the decay rate of elec- 
tronic conductance has been predicted to fall as a result 
of the increasing weight of long-lived, narrow-linewidth 
states • The leading correction to the diffusion predic- 
tion for the electron survival probability P s (t) was calcu- 
lated by Mirlin using the supersymmetry approach 
[iH to be, - hxP s (t) = (*/ti)(1 - t/2iT 2 gTi), where g is 
the dimensionless conductance, which can be expressed 
as g = bv / Av. 

An ideal way to investigate the applicability of the 
diffusion model to mesoscopic systems is to consider 
pulsed electromagnetic transmission. Previous studies 
Hi E H 03 have found exponential decay for t > n, 
as predicted by diffusion theory. However, measurements 
of optical transmission [J^j indicate that the pulse rises 
earlier than predicted by diffusion theory. Even more 
puzzling is the finding by Kop et al. 15] of an increase 
in the inferred value of the diffusion coefficient with in- 
creasing L. 

In this Letter, we report a dramatic breakdown of dif- 
fusion in microwave measurements in nominally diffusive 
random samples for which bv > Av. We find that the 
decay rate of pulsed transmission falls nearly linearly in 
time, as predicted by Mirlin 0. These results are in- 
terpreted in terms of the distribution of decay rates of 
quasimodes of the sample, which is found by taking the 
inverse Laplace transform of the decaying signal. The 
distribution of the modal decay rates is nearly gaussian, 
with an average value that scales as L~ 2 05 , which is close 
to the inverse square scaling for the decay rate in the dif- 
fusion model. The width falls as L -181 , which is faster 
than predicted by Ref. 10]. 

Spectra of the in- and out-of-phase components of the 
steady-state transmitted microwave field are measured 
in low-density collections of dielectric spheres using a 
Hewlett-Packard 8772C vector network analyzer. These 
spectra are multiplied by a gaussian envelope of width 
Af centered at f c and then Fourier-transformed to give 
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the response to a gaussian pulse in the time domain. Cir- 
cular horns are positioned 30 cm in front of and behind 
the sample. Linearly polarized microwave radiation is 
launched from one horn and the cross-polarized compo- 
nent of the transmitted field is detected with the other 
to eliminate the ballistic component of radiation. The 
sample is composed of alumina spheres of diameter 0.95 
cm and refractive index 3.14, contained within a 7.3-cm- 
diameter copper tube at an alumina filling fraction of 
0.068 0. This low density is produced by embedding 
the alumina spheres in Styrofoam spheres of diameter 1.9 
cm and refractive index 1.04. Measurements for random 
ensembles are obtained by momentarily rotating the tube 
about its axis to create new random configurations before 
each spectrum is taken. In this way, measurements are 
carried out in ensembles of 10,000 sample realizations at 
lengths of 61, 90, and 183 cm. In addition, measurements 
are made in an ensemble of 2,300 realizations of a more 
strongly absorbing sample of 90 cm length (L = 90* cm) 
produced by covering 40% of the inside surface of the 
tube with a thin strip of titanium foil laid from end to 
end. The measurements are made within the frequency 
interval 14.7-15.7 GHz for L = 61, 90 and 90* cm and 
15.0-15.4 GHz for L = 183 cm. The frequency inter- 
vals are chosen to be far from sphere resonances |lfij . so 
that the dynamics of transmission is uniform over the 
frequency range and the sample is far from the localiza- 
tion threshold. The closeness to localization, even in the 
presence of absorption, is indicated by the variance of 
the steady-state transmitted intensity normalized to its 
ensemble average value, var {1/(1}) UM- In the absence 
of absorption, var{I/(I)) ~ 1 + 4/3g [H [3. At the 
localization threshold, var(I/(I)) ~ 7/3. The values of 
var(I/(I)) in the samples of L = 61, 90*, 90, and 183 cm 
are 1.18, 1.25, 1.26, and 1.50, respectively. 

The envelope of the temporal response to a gaussian 
pulse peaked at t = is squared to give the transmitted 
intensity I(t) for a particular sample realization. The av- 
erage transmitted intensity (I(t)) is found by averaging 
over the ensemble, and then over the frequency interval 
by shifting f c . The results are shown on a logarithmic 
scale in Fig. la. We find that, when Af > 6v, the tail of 
(/(*)) does not depend on Af. We use Af = 15 MHz for 
L = 61, 90, and 90* cm and Af = 7.5 MHz for L = 183 
cm, so that Af > Sv in all cases. The noise in the fre- 
quency spectra produces a constant background intensity 
in the time domain. Once this background is subtracted, 
a dynamic range of more than 6 orders of magnitude 
is achieved. This makes it possible to study transmis- 
sion on time scales an order of magnitude longer than 
the time of peak transmission, though the longest times 
are still smaller than the inverse level spacing, 1/Av, 
which gives the time required for a photon to explore 
each coherence volume of the sample. The measured de- 
cay rate due to leakage out of the sample and absorption, 
7 = —dln(I(t))/dt, is plotted in Fig. lb. This rate is not 
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FIG. 1: (a) Average transmitted intensity in the alumina sam- 
ples of L = 61 (red), 90* (purple), 90 (blue), and 183 cm 
(green), (b) Temporal derivative of the logarithm of the in- 
tensity gives the rate 7 of the intensity decay due to leakage 
out of the sample and absorption. The thin black curves are 
the decay rates of the diffusion model with D* , zo, and r a 
from Table I. 



constant as predicted by diffusion theory, but falls nearly 
linearly with time. The constant increase in the decay 
rate for the L — 90* sample over that for the L = 90 
sample, seen in Fig. lb, indicates that the absorption 
rate l/r a is constant in time. The decrease in 7 with 
time is thus attributable solely to propagation out of the 
sample, which proceeds at a rate, w 2 D(t) / (L+2zq) 2 . The 
absorption rate is found from a fit of the diffusion model 
[I2I IliL llij to measurements of (I(t)) up to the time 
at which 95% of the full pulse energy has been trans- 
mitted (Fig. 2). The fit is obtained by taking D(t) to 
be a constant, D* , and by minimizing the parameter 
X 2 = (0 _1 E[(J(*)>i - HU)?, where (/(«)), are the 
values of the measured intensity, a is the uncertainty in 
(I(t))i, averaged over the time of the fit, and -T(ij) are 
the values of the model intensity Q| calculated at £j. 
The parameters D* , r a , and zq obtained from the fit are 
listed in Table I and the corresponding decay rates are 
shown by thin solid lines in Fig. lb. For L — 61 cm, 
X 2 at the minimum depends only weakly on r a since the 
temporal range used in the fit is smaller than the r a . For 
this reason, we use the value r Q = 97 ns obtained in the 
fit to the data for L — 183 cm for this length. 

The diffusion coefficient D* obtained from the fit is 
found to decrease slightly with increasing L. The equal- 
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FIG. 2: Fit of the diffusion model to 95% of the full pulse 
energy (circles) transmitted in the alumina samples of L = 61 
(b), 90* (c), 90 (d), and 183 cm (e). The solid curves show 
the fit. (a) is the incident pulse of Af = 15 MHz used to 
obtain (b), (c), and (d); a pulse of Af = 7.5 MHz was used 
to obtain (e). All the curves are normalized. 



ity of the values of D* and zq for the samples with L = 90 
and 90* cm, seen in Table I, indicates that these param- 
eters are not sensitive to the absorption rate. Moreover, 
when 1/T a is subtracted from 7 to give the decay rates 
without absorption in Fig. 3a, the same time-dependent 
decay is found for both samples with length of 90 cm. 

The leakage rates in Fig. 3a give the "time-dependent 
diffusion coefficient", D(t) = (7 - l/r Q )(£ + 2zq) 2 /ti 2 . 
The values of D(t) found from the measurements may 
be compared to those from the theory of Ref. for the 
survival probability, which yields D(t) = D(l — t/ir 2 gTi + 
. . .). This is done by plotting D(t) as a function of the 
dimensionless time t' — tjgr\ in Fig. 3b. The time gr% 
is proportional to l/Av and all the data are predicted 
to fall on a single curve. Although the data for sam- 
ples with different values of L at a given time appear to 
coincide within the noise, there is a clearly discernable 
decrease in the slope of D(t') as L increases. As a result, 
the curves do not extrapolate to a constant bare diffu- 
sion coefficient at t = 0. When D(t) is plotted instead 
versus the dimensionless time t" = t/^/gri in Fig. 3c, the 
slope of D(t") appears to be the same for all values of 



TABLE I: Values of the diffusion coefficient D* , absorption 
time r a , and extrapolation length zo, obtained from fitting 
Eq. (1) of [l^l to the short-time transmitted intensity in Fig. 2. 



L (cm) 


D* (cm 2 /ns) 


Ta (ns) 


z (cm) 


61 


39.4 ± 0.3 


[97] 


9.6 ± 0.3 


90 


37.9 ± 0.3 


104 ± 7 


9.8 ± 0.6 


90* 


37.4 ± 0.4 


46 ± 2 


8.7 ± 0.8 


183 


37.0 ± 0.8 


97 ± 4 


12.1 ± 2.5 



800 
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FIG. 3: (a) Leakage rates in alumina samples with L = 61 
(red), 90* (purple), 90 (blue), and 183 cm (green). The black 
curves are the best fit to the data by a polynomial of power 
2. (b) "Time-dependent diffusion coefficients", D(t) = (7 — 
1/t o )(L+2z ) 2 /V 2 , plotted versus t' = t/gn. (c) D(t) plotted 
versus t" = t/s/gri. 



L, though the curves do not overlap. A strong deviation 
from exponential decay at t" 1 has also been found in 
numerical simulations |20| . 

The changing slope of D{t) with L is associated with 
the scaling of the width of the distribution P(a) of de- 
cay rates of quasi-normal modes of the sample. These 
are hypothesized to form a complete set 0], even when 
8u > Alt. Since the time evolution is given by the su- 
perposition of these modes 0] , the average transmission 
can be expressed as (I(t)) oc J °° P(a) exp(—at)da, pro- 
vided the mode coupling is, on the average, independent 
of a. To find P{a), we use an approximate Laplace in- 
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FIG. 4: Distribution of the modal decay rates, P(a), in the 
alumina samples of L = 61 (red), 90 (blue), and 183 cm 
(green). Inset: scaling of the average (a) (circles) and of 
the width a a (squares) of P{a), and the corresponding best 
fit to ~ (L + 2zo) /3 (solid) are shown on the log- log scale. The 
best fits are obtained with exponents of —2.05 and —1.81 for 
(a) and Set, respectively. 



version algorithm based on the Weeks method [2l|. A 
polynomial of power 2 is fit to the decay rates in Fig. 3a. 
These fits are then used to compute curves (I(t)), which 
are inverted to obtain the distributions P(a) shown in 
Fig. 4. Note that a linear decrease in ^(t) with time as 
7(i) = a — bt would correspond to a gaussian distribution 
of P(a) with (a) = a and var(a) — a 2 a — b. Since the 
decay of j(t) is nearly linear, P(a) is nearly gaussian. 

The scaling of the average decay rate due to leak- 
age out of the sample, seen in Fig. 4, is given by 
(a) oc (L + 2zo)~ 2 ' 05 . This is close to the inverse 
square scaling of the diffusion model and suggests that 
the dynamics observed is characteristic of extended waves 
and is not associated with the approach to localization 
with increasing L. The width of the distribution scales 
as <j a cx (L + 2z )~ 1 ' 81 . This is close to the scaling, 
a a cx Vb cx \/ g x l A Ti oc L -1 ' 75 , and differs from the scal- 
ing predicted by Ref. ^(J> a a cx \/b cx \j ^fgr\ oc 

The wide distribution of the modal decay rates in thin 
samples may be the source of the sharp spectral peaks 
observed in amplifying random media [22J, |23j . In these 
samples, lasing would occur in the longest-lived modes, 
which have the lowest critical gain [24| . 

In conclusion, we have found nonexponential decay of 
pulsed transmission through disordered media in which 
the level width exceeds the spacing between levels, even 
at long times and in thick samples. This departure from 
diffusion theory is interpreted in terms of the decay rate 
statistics of electromagnetic quasi-normal modes. The 
statistics of these modes is fundamental to understanding 
the static and dynamic behavior of waves in both passive 
and active random media. 
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